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Abstract 

Niche theory in its various forms is based on those environmental factors that permit 
species persistence, but less work has focused on defining the extent, or size, of a 
species' environment: the area that explains a species' presence at a point in space. 
We proposed that this habitat extent is identifiable from a characteristic scale of 
habitat selection, the spatial scale at which habitat best explains species' occurrence. 
We hypothesized that this scale is predicted by body size. We tested this hypothesis 
on 12 sympatric terrestrial mammal species in the Canadian Rocky Mountains. For 
each species, habitat models varied across the 20 spatial scales tested. For six species, 
we found a characteristic scale; this scale was explained by species' body mass in 
a quadratic relationship. Habitat measured at large scales best-predicted habitat 
selection in both large and small species, and small scales predict habitat extent 
in medium-sized species. The relationship between body size and habitat selection 
scale implies evolutionary adaptation to landscape heterogeneity as the driver of 
scale-dependent habitat selection. 



Received: 23 May 201 1 ; Revised: 31 August 
201 1; Accepted: 02 September 201 1 



doi: 10.1002/ece3.45 



Introduction 

Niche theory in its various forms (Chase and Leibold 2003) is 
based on those factors that permit species persistence. How- 
ever, much less work has focused on defining the extent, or 
size, of a species' environment: the area, or spatial scale, that 
predicts a species' presence at a point in space. This scale 
describes the area that integrates an individual's lifetime in- 
teractions with competitors, predators, parasites, and mates 
(e.g., Wiens et al. 1993), and so serves as a spatial basis for 
studies of these ecological mechanisms. 

Attempts to quantify this scale inevitably leads to the "scale 
problem" in ecology, wherein the selection of a scale of exper- 
imentation can radically affect the outcome of an experiment, 
a trait that tends to prevent reliable inference across systems 



(Levin 1992). Fundamental to scale theory are the two di- 
mensions of scale: grain and extent. In ecological processes, 
grain is the limit of an organism's perception of heterogene- 
ity, and extent is the maximum area affecting the organism. 
In analysis, grain refers to the lower limit of resolution, or 
"pixel size"; it is the smallest distinct unit of space (or time) 
being measured. Extent defines the size of the experimental 
unit. Kotliar and Wiens (1990) and Levin (1992) maintain 
that matching the grain and extent of an experiment to the 
grain and extent of the process it measures must be the first 
step in ecological inquiry. With the development of con- 
ceptual frameworks to solve the scale problem, such as scale 
domains (Wiens 1989) and ecological hierarchy theory (Allen 
and Hoekstra 1992), the concept of scale in ecology has de- 
veloped from a pervasive cautionary tale (Levin 1992) to an 
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increasingly effective scientific tool ( Schneider 2001). Exami- 
nation of ecological processes across scales has revealed much 
about these processes, from allometric relationships (Peters 
1983) to biodiversity patterns (Storch et al. 2007). 

For example, multiple-scale analysis of density-dependent 
habitat selection processes (Morris 1987; Morris 2003) 
and species-habitat relationship patterns (e.g., Fisher et al. 
2005; Knegt et al. 2010b) show that habitat selection varies 
markedly across spatial scales. How and why it varies is 
hotly debated, due partly to inconsistencies in experimen- 
tal design, data collection, and analysis among studies and 
species (Schooley 2006; Wheatley and Johnson 2009). It is 
therefore not surprising that recent attempts to derive gen- 
eralities have highlighted the many controversies inherent in 
multiscale habitat selection (Bowyer and Kie 2006; Mayor et 
al. 2009a; Wheatley and Johnson 2009; Laca et al. 2010). It is 
unknown whether scale-dependent habitat selection is driven 
primarily by population-level responses (Morris 2003), indi- 
vidual responses that are a function of home-range size (e.g., 
McLoughlin et al. 2004), or both. 

At the heart of these and other major conceptual gaps lays 
an unanswered question: are several spatial scales required 
to explain a species occurrence at a point in space, or is 
there a single best scale? Several lines of evidence suggest that 
ecological processes operate simultaneously at different do- 
mains of scale; this basic prediction from hierarchy theory 
is the basis for hierarchical habitat selection (Johnson 1980), 
wherein different processes, such as foraging and home-range 
selection, occur at different scales. Notably, as Levin (1992) 
suggested, not all scales are equal. It is possible that within a 
species, one scale may better describe habitat selection than 
other scales. Holland et al. (2004) termed this the character- 
istic selection scale. A characteristic selection scale has been 
demonstrated for fly parasitoids in boreal forests (Roland 
and Taylor 1997), forest beetles in mixedwood forests (Hol- 
land et al. 2004), grizzly bears ( Ursus arctos horribilis; Nams 
et al. 2006), and marten (Martes americana; Mowat 2006). 
Though the existence of a characteristic scale of selection 
within species has been suggested by several empirical stud- 
ies, it remains contentious (Holland et al. 2004; Bowyer and 
Kie 2006; Schooley 2006). 

If characteristic spatial scales of habitat selection exist, 
these may arise from a species' interaction with landscape 
patchiness and resource dispersion. The spatial structure of 
the landscape has been shown to affect the scale of habi- 
tat selection for caribou (Mayor et al. 2007, 2009b, Schaefer 
and Mayor 2007) and elephants (Knegt et al. 2010a, 2010b). 
Species interact with landscape spatial structure through 
combined processes of foraging, dispersal, and migration. 
Since these processes scale with body mass (Peters 1983), we 
hypothesize that the characteristic scale of a species' habitat 
selection increases linearly as a function of its body mass. 
This hypothesis can be tested by examining habitat selec- 
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tion across scales of different-sized species, occurring in the 
same landscape, at the same time. We test this hypothesis 
by simultaneously surveying 12 mammal species sympatri- 
cally distributed in a mountain landscape, and modeling the 
habitat selection of each species across 20 spatial scales. We 
test whether a characteristic habitat selection scale exists for 
these mammal species, and whether this scale is predicted by 
species' body mass. 

Methods 

Study area 

Terrestrial mammals were sampled in the Foothills and Rocky 
Mountains of west-central Alberta, Canada, over a 6,400-km 2 
area. The entire Rocky Mountain portion of our study site lies 
within the Willmore Wilderness Park, a 4,600-km 2 conser- 
vation area protected from forest harvesting, mining, seismic 
exploration, and roads. Topography is rugged; alpine areas 
are characterized by mountain meadows of herbs and shrubs. 
Subalpine slopes are forested by Picea engelmanni and Abies 
lasiocarpa. The Foothills form the eastern border of the Rocky 
Mountains, with moderately rugged topography and lower 
elevation than the Rockies. Forests are most commonly Pinus 
contorta with P. glauca. Forest harvesting and energy devel- 
opment are common. 

Study design 

Baited survey sites were deployed in December and sampled 
monthly until March, a period during which food is scarce 
and bait attracts mammals. Differences in logistical require- 
ments in the Rockies versus Foothills necessitated different 
sampling designs. In the Rockies, we used a systematic sam- 
pling design constrained by helicopter access and avalanche 
risk. Sixty-six sites were placed an average of 5,727 m apart 
(SD = 1,574 m); 30 were sampled in 2006-2007 and the re- 
mainder in 2007-2008. In the Foothills, we deployed sites in 
2005-2006 along a ca. 415-km transect, using a constrained 
systematic design that followed access and omitted wetlands, 
areas adjacent (<1 km) to main highways, and current in- 
dustrial activity. Fifty-four sites were deployed, an average of 
4,335 m apart (SD = 5,218 m), ca. 50 m from access roads. 

Sampling species occurrence 

We used a combination of noninvasive genetic tagging (NGT) 
via hair sampling and camera trapping to survey mammal 
occurrence. NGT samples were collected using Gaucho® 
barbed wire (Bekaert, Brussels, Belgium) wrapped around 
a tree rebaited monthly with a whole beaver (ca. 10 kg). 
Species were identified from hair DNA (Wildlife Genetics 
International, Nelson, British Columbia, Canada). DNA was 
extracted from hairs using QIAGEN®'s DNEasy™ Tissue 
Kits (QIAGEN, Hilden, Germany) and analyzed to identify 
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species using sequence-based analysis of the 16S rRNA gene 
of mitochondrial DNA (mtDNA; sensu Johnson and O'Brien 
1997) that was then compared against a DNA reference library 
of all known mammal species in the region. We summed 
presences across 3 months to yield a 0-3 count of species 
occurrences at each site. NGT captured only mustelids; so in 
the Rockies, we also used infra-red-triggered digital camera 
traps (O'Connell et al. 2010). We deployed Reconyx RM30 
(2006-20087) or PM30 (2007-2008) infrared-triggered dig- 
ital cameras (Reconyx, Holmen, WI, USA) at all 66 sampling 
sites. Camera traps were effective at detecting more mammal 
species that visited a site. We summed each species' presence 
across 3 months to yield a 0-3 count of species occurrences 
at each site. 

Habitat analysis 

We used a LandSat thematic-mapped GIS land cover dataset 
incorporating a digital elevation model classified using a 
habitat-identification algorithm (McDermid et al. 2009). 
This dataset yielded 16 potential habitat types. We retained 
seven that were hypothesized to be of ecological importance 
to the study species, and which occurred sufficiently often 
in the study landscape to allow modeling: closed conifer 
forest, moderate conifer forest, open conifer forest, mixed- 
wood forest, open wetland, upland shrubs, and upland herba- 
ceous habitats (see McDermid et al. 2009 for descriptions). 
The remaining nine variables occurred very infrequently 
in these landscapes and were not expected to be impor- 
tant predictors of occurrence. We used ArcGis 9.3 (Envi- 
ronmental Systems Research Institute Inc., Redlands, CA, 
USA) Spatial Analyst to geo-reference each sampling site on 
the habitat GIS layer. Using spatial analysis routines (writ- 
ten based on Arc- View v3.x Spatial Analyst) and the Re- 
gional Analysis function of Patch Analyst, an extension to 
ArcView (http://flash.lakeheadu.ca/wrrempel/patch), we cal- 
culated the percentage of each land-cover habitat type at 
20 spatial scales, consisting of a circle with radii of 250- 
5,000 m, in 250-m increments. Thus, holding grain constant, 
we analyzed 20 different extents, which we refer to as spatial 
scales. 

Statistical analysis 

We modeled species occurrence against available habitat sur- 
rounding each sampling point using generalized linear mod- 
els (Poisson errors, log link) in R ver. 2.10.1 (R Develop- 
ment Core Team, 2010). Mustelids were analyzed using the 
combined Foothills and Rockies dataset; other species were 
analyzed using only the latter. Each candidate set of mod- 
els consisted of combinations of the seven habitat variables 
measured at a single spatial scale. For each candidate set of 
models - one set for each scale - we used the step-AIC func- 
tion (MASS package for R software) to identify the best-fit 
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model, defined as the model with the lowest (Akaike infor- 
mation criterion) AIC score (Burnham and Anderson 2002). 
AIC weights were calculated for each best-fit model, and plot- 
ted against the spatial scale at which that model was derived. 

Recognizing that each spatial scale we chose may only ap- 
proximate the interval in which the possibly "true" character- 
istic scale lies, we defined a characteristic scale as a single scale, 
or range of spatial scales, at which the species-habitat model 
is supported by >50% of the weight of evidence as defined 
by summed model AIC weights. We repeated this method 
for each species to determine whether a characteristic scale 
of selection could be identified. Where a characteristic scale 
was identified, we tested whether this scale was significantly 
explained by body mass and home-range size values (from 
Holling 1992) using multiple linear regression models. 

To quantify patterns in spatial structure of our landscapes 
across scales, we calculated the variance of each habitat metric 
at each spatial scale and plotted these across scales (sensu 
Wheatley2010). 

Results 

We obtained repeat detections of Glaucomys sabrinus (fly- 
ing squirrel), Tamisciurus hudsonicus (red squirrel), Lepus 
americanus (snowshoe hare), Mustela erminea (ermine), M. 
americana (marten), M. pennanti (fisher), Gulogulo (wolver- 
ine), Vulpes vulpes (red fox), Canis latrans (coyote), C. lupus 
(wolf), Lynx canadensis (lynx), and Puma concolor (cougar). 
Mustelids were detected in both the Foothills and the Will- 
more, whereas all other species were only detected in the 
Willmore, where we had deployed infrared-triggered cam- 
eras. Detection frequency («) differed among species, and 
consequently so did model fit (Table 1). 

Peaks in characteristic scales 

Cougars, wolves, coyotes, and flying squirrels were observed 
at fewer than 20% of the sampling sites. Analysis of model 
AIC weights across scales for these species did not reveal peaks 
for any species but cougar. Among the eight species detected 
at >20% of the sites, we observed peaks in selection scale for 
six. 

Among sciurids, red squirrels showed a definitive peak 
between 1,750 m and 2,250 m radii (Fig. 1), supported by 
62.3% of the weight of evidence. Flying squirrels showed no 
pattern. The only lagomorph — snowshoe hares — showed a 
very weak multimodal pattern. Among mustelids, wolverine 
showed the clearest peak (Fig. 2); 64.6% of the weight of 
evidence suggests habitat measured at the 4,750-5,000 m 
scales best explains occurrence. Fisher also showed a clear 
pattern, as 81.4% of the weight of evidence suggests 500- 
m radius is the best selection scale. For ermines, 61.5% of 
the weight of evidence supported a characteristic scale of 
4,500-5,000 m. Marten showed a much weaker pattern, as 
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Table 1. Generalized linear models (Poisson errors) of mammal occurrence against percent cover of habitat types in the surrounding landscape. 
Habitat was measured at 20 different spatial scales; the best-supported model (highest AIC weight) among 20 models is shown for each species. 



Species common 






Characteristic 


Null 




Residual 




name 


n* 


Mass (g)** 


scale 


deviance 


df 


deviance 


df 


Cougar 


8 


66,508 


250 


43.29 


65 


24.72 


62 


Wolf 


8 


43,205 


3,000 


38.64 


65 


18.01 


61 


Flying squirrel 


9 


105 


4,500 


50.28 


65 


38.10 


63 




1 1 


1 4 061 




55 54 


65 


34 61 


63 


Lynx 


15 


10,149 


1,250 


66.82 


65 


47.94 


62 


Ermine 


17 


81 


5,000 


97.47 


65 


63.54 


63 


Hare 


19 


1,497 


5,000 


72.54 


65 


58.78 


62 


Fox 


21 


5,193 


250 


77.72 


65 


63.56 


62 


Red squirrel 


39 


191 


2,000 


77.80 


65 


58.31 


63 


Fisher 


48 


3,118 


500 


188.92 


119 


147.88 


115 


Wolverine 


48 


12,303 


5,000 


183.58 


119 


104.94 


116 


Marten 


62 


839 


4,500 


164.75 


119 


118.26 


115 



n* = number of sites at which the species was detected. 
"Masses are rounded from Holling (1992). 



AIC weights did not exceed 0.11. Weights increased from the 
2,500-m scale to the 5,000-m scale, but no distinct peak was 
evident. 

Among canids, red fox habitat selection showed a distinct, 
high peak at the 250-m scale (Fig. 3), supported by 66.3% 
of the weight of evidence. Coyote habitat selection showed a 
less distinct pattern. AIC weights increased from 2,750-m to 
4,250-m scales, then decreased, though AIC weights did not 
exceed 0.14. Wolves showed a similar weak pattern, except 
that AIC weights increased from 1,750 to 3,000 m, and then 
decreased. Both felids showed distinct peaks in habitat selec- 
tion scale (Fig. 4). Over half the weight of evidence (51.6%) 
suggested lynx selected habitat at scales between 1,000 m and 
1,250 m. Habitat models for cougars peaked at 250 m (AIC 
weight = 0.378) and no other models were supported, but low 
detection and poor model fit (Table 1) suggests this model is 
spurious. 

Variance in habitat metrics 

Variance in percent cover of each habitat variable decreased 
with spatial scale (Fig. 5). Variance was highest at the smallest 
scale; in general, the variance decreased rapidly to the 1,000- 
m scale, and then reached an asymptote. This pattern was 
similar for all seven habitat variables used in the analysis. 
Maxima or minima in variance did not match characteristic 
habitat selection scales across species. 

Relationship between body mass and 
habitat selection scale 

For the six species with a defined characteristic scale of 
selection — red squirrel, ermine, wolverine, fisher, red fox, 
and lynx — home-range size (from Holling 1992) was not a 



significant predictor of characteristic selection scale (linear 
model; P = 0.216; F = 2.16; df = 1 and 4; adjusted R 2 = 
0.35). Contrary to our hypothesis, body mass was also not a 
significant predictor in the linear model (P = 0.853; F = 0.04; 
df = 1 and 4; adjusted R 2 = 0.0097). However, the pattern 
of the plotted data led us to test whether the relationship be- 
tween body mass and selection scale may be logarithmic and 
quadratic, so we modeled log(selection scale) against mass 
and mass-squared. This analysis suggested that a quadratic 
function of body mass is a significant predictor of selection 
scale. Selection scale decreased with increasing body mass, 
and then increased again (Fig. 6). Model terms are statisti- 
cally significant (mass P = 0.0065; mass 2 P = 0.0056; model 
P = 0.0128; F = 25.97; df = 2 and 3). The model explains 
a high percentage of variation in the data (multiple R 2 = 
0.95). 

Discussion 

This first explicit test for characteristic scales of selection 
across species demonstrates a link between body size and 
spatial scales of habitat selection. Several key points emerge. 
First, a characteristic habitat selection scale was definable 
for half of the 12 mammal species we studied. The marked 
peaks in model support at distinct scale domains for these 
species strongly suggest the existence of a characteristic scale 
of habitat selection, similar to Holland et al. (2004), Nams 
et al. (2006), and Mowat (2006). 

The existence of a characteristic spatial scale of selection 
has some key theoretical implications. Our results support 
Wiens' (1989) and Levin's (1992) hypothesis that ecologi- 
cal processes operate at characteristic scale domains. Though 
widely accepted, there have been few empirical tests of this 
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Figure 1. Support for models (AIC weights) of sciurid and lagomorph occurrence against habitat measured across 20 spatial scales (circles of 250- 
5,000 m radii) around each sampling point. 



hypothesis across a broad range of spatial scales (but see Hol- 
land et al. 2004; Nams et al. 2006). Our results also support 
their contention that not all scales serve equally well as a 
basis for examining ecological processes; some scales are bet- 
ter than others, since ecological processes manifest at some 
scales but not others. For example, models relating red squir- 
rel (a conifer specialist) occurrence to open conifer forest at 
the 250-m scale yielded a nonsignificant relationship (P = 
0.0963, null deviance = 77.8, residual deviance 74.6, df = 65). 
This same model at the 2000-m scale yielded a significant re- 
lationship (P = 0.000685, residual deviance = 62.3). If we 
test for habitat selection at the wrong scale, we fail to detect 
existing relationships that become observable at appropri- 
ate scales; this remains an ongoing problem, despite greater 
acceptance of scale dependency in habitat selection (reviewed 
in Wheatley and Johnson 2009). This is arguably true of all 



studies of scale-dependent ecological processes: study scales 
must be empirically modeled and matched to process scales, 
or risk flawed inferences and conclusions. 

We found that characteristic habitat selection scales vary 
among species, suggesting that some scale-dependent ecolog- 
ical mechanism is driving habitat selection differently for each 
species. The importance of landscape structure for species has 
also been indicated by Roland and Taylor (1997), Nams et al. 
(2006), Mowat (2006), Schaefer and Mayor (2007), Mayor 
et al. (2009b), and Knegt et al. (2010a), but the mechanism is 
still debated. If we had observed that the same habitat vari- 
ables were significant for small and large mammal species, 
we might suspect that a shared response to landscape pat- 
tern was driving the observed curvilinear relationship. This 
was not the case; different combinations of variables emerged 
as significant for each species. Alternatively, the variance in 
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Figure 2. Support for models (AIC weights) of mustelid occurrence against habitat measured across 20 spatial scales (circles of 250-5,000 m radii) 
around each sampling point. 



habitat metrics may be expected to predict the characteris- 
tic scale of habitat selection (e.g., Mayor et al. 2007, 2009b; 
Schaefer and Mayor 2007; Wheatley 2010; Knegt et al. 20 1 0a) . 
However, variance of the habitat metrics decreased from the 
250-m scale to the 5,000-m scale. Scales with maximum or 
minimum variance are not the scales at which habitat best 
predicted species occurrence. Though further investigation 
with different landscape metrics may add insight, our results 
do not support the hypothesis that scale-dependent changes 
in habitat and landscape structure match species' character- 
istic selection scales. 

Instead, our analysis suggests that characteristic scales of 
habitat selection are related to body size. Body size has been 
previously implicated as a driver of spatial scales of herbivore 
foraging (Wilmshurst et al. 2000; Laca et al. 2010) due to 
constraints on rates of food intake imposed by body size 
(Peters 1983), and has a well-known relationship with home- 



range size (Haskell et al. 2002). However, all but three of 
our study species are carnivores, which do not forage directly 
on the vegetation within the habitat types we quantified. 
This suggests the scale-dependent response to habitat extends 
beyond individuals' foraging patterns within and between 
food patches. Rather, species are likely responding to several 
concurrent ecological processes structured by heterogeneous 
landscapes {sensu Wiens et al. 1993). We suggest that the 
spatial scale of habitat selection is the result of a species' 
interaction with the structure of its landscape, mediated by 
its body size. 

Our conclusions are consistent with Holling's (1992) 
textural-discontinuity hypothesis. Holling (1992) suggested 
that body-size gaps represent scale domains wherein 
resource-patch size is too small to provide enough resources 
to support a body mass x, but patch dispersion is too great for 
a stride length f(x) to efficiently use multiple patches. Holling 
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Figure 3. Support for models (AIC weights) of canid occurrence against habitat measured across 20 spatial scales (circles of 250-5,000 m radii) 
around each sampling point. 
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Figure 4. Support for models (AIC weights) of felid occurrence against habitat measured across 20 spatial scales (circles of 250-5,000 m radii) around 
each sampling point. 
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Figure 5. Variance in habitat metrics (percent cover of patch type) across spatial scales. Variance is highest at small scales as they tend to truncate 
average patch size, then decreases as scale increases. 



also suggested that within a scale domain, a species either 
had to be small enough to extract energy from resources at- 
tainable within the limits of its vagility, or large enough to 
use multiple patches. Species evolve body masses aggregating 
about a mean that is adapted to the characteristic resource 



patchiness within each scale domain (Holling 1992). Thus, 
the probability of a species' occurrence in space is explained 
by a spatial extent related to its body size. 

If body-size distributions arose from evolutionary con- 
straints imposed by discontinuities in spatial landscape 
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Figure 6. Characteristic scale of habitat 
selection (determined by AIC weight, see 
Figs. 1 and 2), log-transformed and modeled 
against body mass of six mammal species for 
which a characteristic scale was detectable. 
Habitat quantified at large scales best predicts 
both small and large mammal occurrence, 
whereas habitat quantified at small scales best 
predicts occurrence of intermediate-sized 
mammals. 
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patterns, then body size may predict species' scale-dependent 
response to landscape structure. We expected such a rela- 
tionship between spatial scale of selection and body size to 
be linear, reflecting (for example) the ecological processes of 
landscape complementation and supplementation (Dunning 
et al. 1992). The quadratic relationship between body mass 
and selection scale that we observed was unexpected. We hy- 
pothesize that the quadratic relationship reflects individual- 
level responses to landscape structure for large mammals, 
and population-level responses for small mammals. By way 
of analogy, in large mammals we have detected a behav- 
ioral response, whereas in small mammals we have detected a 
numerical response. Within the range of spatial extents that 
we measured, mid-size and large mammals select denning 
and shelter sites, foraging patches, and home ranges based 
on resource availability (sensu Johnson 1980). The spatial 
scale at which this selection occurs is mediated by species' 
body size. So, mid-size mammals select habitat that reflects 
resource requirements at smaller scales than do larger mam- 
mals. In contrast, for small mammals we hypothesize that 
landscape heterogeneity affects population processes such as 
dispersal, immigration, and emigration, rather than forag- 
ing (Pulliam 1988). Patch size, proximity, and heterogeneity 
drive population dynamics and determine probabilities of 
local population persistence at large scales. Individual-level 
resource selection likely occurs at a smaller grain (higher 
resolutions) than we measured. 



The dichotomy in the ecological processes creating the ob- 
served patterns in scale-dependent habitat selection might 
result in the quadratic relationship we measured. Large 
scales predict small mammal occurrence through popula- 
tion processes; and the scale of habitat selection, through 
individual-level resource selection, competition, and pre- 
dation, increases with increasing body size from mid-size 
through larger mammals. Some support for this hypothesis 
may already exist. Recently Brown (2007) found that body 
size mediated birds' response to forest habitat patch size and 
isolation. Very small and very large birds responded similarly 
to landscape structure, whereas large and mid-sized birds re- 
sponded differently. Brown (2007) attributed this pattern to 
different dispersal and space-use abilities among size classes, 
analogous to the limitations imposed by body size and stride 
length in nonvolant mammals. 

We did not find a characteristic scale of selection for some 
species, for which we suggest the following potential expla- 
nations. (1) There were insufficient observations to detect 
a characteristic scale. This is most likely to apply to flying 
squirrels, coyote, wolf, and cougar. (2) While we have used a 
range of spatial scales to look at these relationships, we recog- 
nize that these scales could be too large for some species (e.g., 
flying squirrel) and not large enough for others (e.g., coyote, 
wolf, and cougar). (3) Some animals have multiple scales of 
selection. We contend that failing to detect an existing charac- 
teristic scale for species with low observation number (Type 2 
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error) is more likely than falsely detecting a nonexistent scale 
for six species (Type 1 error). The lack of a peak for marten 
and hares, however, requires further investigation. It is pos- 
sible that a species might be selecting for different habitat 
variables at different scales. In fact, we would expect this to 
be the case (sensu Johnson 1980), as habitats are distributed 
differently within the study landscape, and their value likely 
varies with their availability and distribution. We could not 
include all variables at all different scales in a single model; 
with n = 120 (wolverines, fisher, marten) or n = 66 (all other 
species), the number of potential variables exceeds n, which 
Burnham and Anderson (2002) strongly recommend against. 
Instead, we modeled variables at a subset of scales (250; 500; 
1,000; 2,000; and 4,000 m) to examine whether different 
variables were significant at different scales (Appendix 1). 
Most models failed to achieve convergence, but those that did 
(marten, fisher, wolverine, red squirrel) produced different 
significant variables at different scales as anticipated, and did 
not necessarily predict species' occurrence better than single- 
scale models (Appendix 1 ) . We might expect scale-dependent 
selection among different variables to reduce our power to 
find a pattern across scales, but this was not the case. The fact 
we did find a characteristic scale of selection for some species, 
despite selection of different habitat variables among scales, 
shows the relationship is robust. 

Our results illustrate that ecological studies examining the 
many environmental variables driving habitat selection must 
be extremely cautious about spatial scale and the area that 
they define as habitat. Our results also suggest that habi- 
tat selection may have evolved in conjunction with body 
size; likewise, body size may have evolved to match the 
distribution of resources on the landscape {sensu Holling 
1992). This relationship provides ecologists a basis for iden- 
tifying the range of scales at which to conduct studies 
of habitat selection and a conceptual foundation for ex- 
amination of the mechanisms driving species' distribution 
patterns. 
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Appendix 1 

Comparison of single-scale and multi-scale habitat selection models for mammal species in mountain regions of Alberta, Canada. The single-scale 
model is that with the lowest AIC score from a candidate set of single-scale models, each with habitat quantified at a spatial scale ranging from 250-m 
to 5000-m radii, at 250-m intervals. The multi-scale model is that with the lowest AIC score from a stepwise selection of a global model containing 
all habitat variables measured at 5 scales - 250 m, 500 m, 1 000 m, 2000 m, and 5000 m. 

Single-scale models Multi-scale models 



Habitat Variables+ Null/Residual AIC Habitat Null/Residual AIC 

Species (Scale) Deviance* Score Variables+ Deviance* Score 



wolverine (5000-m) 


183.6(119) 281.86 HERB 250 


183.6(119) 286.47 


DCON 


105.0 (1 16) SHRUB 500 


95.6 (109) 


HERB 


HERB 500 




MFOR 


MCON 500 






OWET 500 






DCON 500 






MCON 1000 






MCON 2000 






MCON 4000 






MFOR 4000 




fisher (500-m) 


188.9(119) 277.73 MCON 250 


188.9(119) 268.00 


DCON 


147.9 (1 1 5) DCON 250 


120.1 (106) 


SHRUB 


SHRUB 500 




MCON 


OWET 500 




OWET 


DCON 500 






SHRUB 1000 






MCON 1000 






MFOR 2000 






MCON 2000 






MCON_4000 






OCON_4000 






DCON_4000 






MFOR_4000 




marten (4500-m) 


164.7 (1 19) 357.53 SHRUB_250 


164.7(119) 361.36 


OCON 


118.3(115) OWET_250 


108.1 (108) 


HERB 


MCON_500 




OWET 


SHRUB_500 




MFOR 


DCONJ000 






MCONJ000 






OCON_2000 






MCON_2000 






SHRUB_2000 






MFOR_2000 






HERB_4000 




red squirrel (2000-m) 


77.8 (65) 155.01 DCON_2000 


77.8 (65) 152.52 


DCON 


58.3 (63) OCON_2000 


53.8 (62) 


OCON 


DCON_4000 




•(degrees of freedom) 






+DCON = dense conifer forest; MCON = 


moderate conifer forest; OCON = open conifer forest; MFOR = mixed deciduous / conifer forest; OWET = 


open wetland; SHRIB = shrub-dominated cover; HERB = herbaceous-dominated cover. Descriptions in McDermid ef al. (2009). 
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